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A measurement for studying the mass dependence of the dilepton interferometry in relativistic 
heavy-ion collision experiments as a tool to characterize the quark-gluon phase is proposed. In 
calculations involving dileptons, we show that the mass dependence of radii extracted from the 
virtual photon (dilepton) interferometry provide access to the development of collective flow with 
time. It is argued that the non-monotonic variation of HBT radii with invariant mass of the 
lepton pairs signals the formation of quark gluon plasma in heavy ion collisions. Our proposal 
of experimentally measuring the ratio, R ou t /Rside for dileptons can be used to estimate the average 
life times of the partonic as well as the hadronic phases. 
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I. INTRODUCTION 

The primary goal of the nuclear collision experi- 
ments at ultra-relativistic energies is to create and 
study a new state of matter called the Quark Gluon 
Plasma (QGP) 1]. Most of the experimental ob- 
servables for QGP, however, contain contributions 
from both the partonic as well as hadronic phases. 
As a consequence the disentanglement of the signals 
for the partonic phase remains a challenge and has 
been successful to a limited extend [l| . Some of the 
probes, which are produced early in the interactions 
like jets Q and heavy quarks [3j or real photons 
and dileptons (7l-fl0j are considered to be particularly 
useful. 

The transverse momentum (j>t) distribution of 
photons reflect the temperature of the source as their 
productions from a thermal system depend on the 
temperature (T) of the bath through the thermal 
phase space factors of the participants of the reac- 
tion that produces the photon. However, the ther- 
mal phase space factor may be changed by several 
factors - e.g. the transverse kick due to flow re- 
ceived by low pr photons from the low temperature 
hadronic phase will mingle with the high pt pho- 
tons from the partonic phase, making the task of de- 
tecting and characterizing QGP more difficult. For 
dilepton the situation is, however, different because 
in this case we have two kinematic variables - out of 
these two, the px spectra of lepton pairs is affected 
but the pr integrated invariant mass (M) spectra 
is unaltered by the flow. Moreover, it is expected 
that the large M pairs originate from the early time 
and the low M pairs predominantly produced in the 
late time. Therefore, the M distribution can act as 
a chronometer of the heavy ion collisions. This sug- 
gests that a simultaneous measurement of pr and M 
spectra along with a judicious choice of pr and M 
windows will be very useful to characterize the QGP 
and the hadronic phases. Precise measurements of 



lepton pairs in pp collisions at a given collision en- 
ergy is of paramount importance for detecting the 
thermal spectra in heavy ion collisions at the same 
energy. 

Experimental measurements of two-particle inten- 
sity interferometry has been established as a useful 
tool to characterize the space-time evolution of the 
heavy- ion collision [ll|. For the case of dileptons, 
such an interferometry needs to be carried out over 
dilepton pairs, theoretically representing a study of 
the correlations between two virtual photons. Al- 
though, the dilepton production rate is down by a 
factor of a compared to real photon, the analysis 
involving lepton pairs has been successfully used to 
get direct photon yields at RHIC [llj]. In contrast 
to hadrons, two-particle intensity interferometry by 
using lepton pairs, like photons, have almost no 
interaction with the surrounding hadronic medium 
hence can provide information on the history of the 
evolution of the hot matter very efficiently. From 
the experimental point of view, dilepton interferom- 
etry encounters considerable difficulties compared 
to hadron interferometry due to small yield of the 
dileptons from the early hot and dense region of 
the matter and the associated large background pri- 
marily from the electromagnetic decay processes of 
hadrons at freeze-out. However, recent work demon- 
strate that it is still possible to carry out experimen- 
tally such interferometry studies jlSj . With a high 
statistics data already collected at RHIC in the year 
2010 by both STAR and PHENIX collaborations 
having dedicated detectors (Time-Of-Flight [l4[ and 
Hadron Blind Detector [15]) with good acceptance 
for dilepton measurements, also augurs well for the 
dilepton interferometry analysis. 

In this work we present this new proposal for car- 
rying out an experimental measurement of dilepton 
interferometry at RHIC. We establish through a hy- 
drodynamical model based space-time evolution for 
central 0-5% Au+Au collisions at ^Jsnn = 200 GeV 
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the promise of such a dilepton interferometry analy- 
sis will hold out to understand the properties of the 
partonic phase. 

II. DILEPTON INTERFEROMETRY IN 
HEAVY ION COLLISION 

As interferometry of the dilepton pairs actually 
reflect correlations between two virtual photons, the 
analysis then concentrates on computing the Bosc- 
Einstein correlation (BEC) function for two identical 
particles defined as, 

c a (*U)= P2 fc fcl l (i) 

fM*i)iM*a) 
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where K = (fci + k 2 )/2, Afc M = - k 2 ^ = q^, 
x and k are the four co-ordinates for position and 
momentum variables respectively and ui(x, k) is the 
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The inclusion of the spin of the virtual photon will 
reduce the value of C 2 — 1 by 1/3. The correla- 
tion functions can be evaluated for different average 
mass windows, (M) ( = Af e+e -) = (A/i+M 2 )/2. The 
leading order process through which lepton pairs are 
produced in QGP is qq — > [lg]. For the low 

M dilepton production from the hadronic phase the 
decays of the light vector mesons p, uj and 6 have 
been considered including the continuum [114H [I?} • 
Since the continuum part of the vector meson spec- 
tral functions are included in the current work the 
processes like four pions annihilation [l8| are ex- 
cluded to avoid double counting. For the space time 
evolution of the system relativistic hydrodynamical 
model with cylindrical symmetry [19j and boost in- 
variance along the longitudinal direction [2(| has 
been used. The initial temperature (Tj) and proper 
thermalization time (tv) of the system is constrained 
by hadronic multiplicity (dN/dy) as dN/dy ~ Tfn. 
The equation of state (EoS) which controls the rate 
of expansion/cooling has been taken from the lat- 



where fej is the three momentum of the particle i 
and -Pi(fci) and P 2 {k\ 1 k 2 ) represent the one- and 
two- particle inclusive lepton pair spectra respec- 
tively. The correlation function, C 2 has been evalu- 
ated with the following source function: 



Pi(jfc) = / d A x uj(x,k) (2) 



and 



(3) 



source function related to the thermal emission rate 
of virtual photons per unit four volume, R as follows: 
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TABLE I: Values of the various parameters used in the 
relativistic hydrodynamical calculations. 
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tice QCD calculations [2l|. The chemical (T c h) and 
kinetic (T/ Q ) freeze-out temperatures are fixed by 
the particle ratios and the slope of the pr spectra 
of hadron 22]. The values of these parameters are 
displayed in Table HI 

With all these ingredients we evaluate the cor- 
relation function C 2 for different invariant mass 
windows as a function of q S ide and q ou t which are re- 
lated to transverse momentum of individual pair as 
follows: q out = (kj T - fc| T )//(fciT, k 2 r) and q side = 



P 2 (k u k 2 )=P 1 (k 1 )P 1 {k 2 )+ I d 4 x 1 d i x 2 u J (x 1 ,K)u(x 2l K) cos(A^A^) 
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FIG. 1: (color online) Correlation function for dilepton 
pairs as a function of q S ide [(a), for k\r = k-zr = 2 GeV 
and i/)2 = 0] and q ou t [(b), for ipi = ifa = and kyr = 2 
GeV] for three values of (M). The solid lines shows the 
parameterization of C2 using Eq. O 



(2fci T fc 2 T \/l — c °s ; 
f{k 1T ,k 2T ) = V^? 
fcj = (MiTCOshyi,kiTCOs%pi,k iT sin^i,M lT sinhyi) 
and M^t 



- V> 2 ))//(fciT, ^2t) where 
1T + fe|r + 2k 1T k 2 T cos(-0i — V'2), 



< M > 2 



III. RESULTS AND DISCUSSIONS 

Fig. [T] shows the BEC function of dilepton pairs 
for different values of (JW) as a function of qWe and 
<Zout- We have evaluated the C 2 for (M) =0.3, 0.5, 
0.77, 1.2, 1.6 and 2.5 GeV, however, in Fig. Q] we 
display results for only three values of (M) corre- 
sponding to low and high mass which are expected 
to be dominated by radiations (see Fig. [5]) from 
QGP ((M) - 0.3,1.6 GeV) and hadronic phase 
((M) ~ 0.77 GeV) respectively. A clear difference is 
seen in C 2 for different (M) windows when plotted 
as a function of q s %de- The differences are, however, 
small when BEC is studied as a function of q ou t- 

The source dimensions can be obtained by param- 
eterizing the calculated correlation function of the 
(time like) virtual photon with the empirical (Gaus- 
sian) form: 



C 2 



1 + Acxp(-i?^ 



(5) 



where the subscript i stand for out and side and A (= 
1/3 here) represents the degree of chaoticity of the 
source. The deviation of A from 1/3 will indicate the 
presence of non-thermal sources. A representative fit 
to the correlation functions are shown in Fig. Q] (solid 
lines). 

While the radius (R S ide) corresponding to q S ide 
is closely related to the transverse size of the sys- 
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FIG. 2: (color online) Invariant mass distribution of lep- 
ton pairs from quark matter and hadronic matter. 
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FIG. 3: (color online) R 3 ide and R ou t as a function of 
(M). The dashed, dotted and the solid line (with aster- 
isk) indicate the HBT radii for the QGP, hadronic and 
total dilepton contributions from all the phases respec- 
tively. The open circles are obtained by switching off the 
contributions from p and w. 



tern and considerably affected by the collectivity, 
the radius (R ou t) corresponding to q out measures 
both the transverse size and duration of particle 
emission [lj [23[ . The extracted R S id e an d Rout 
for different (M) are shown in Fig. [31 The R S id e 
shows non-monotonic dependence on M. The R S ide 
reduces with (M), reaches its minimum value at 
(M) ~ m p and then increases again at high (M) 
approaching values close to the corresponding R S id e 
for the QGP phase. It can be shown that R S ide ~ 
1/(1 + Ecoiiective/ Ethermai) ■ In the absence of radial 
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flow, R S ide is independent of M. With the radial 
expansion of the system a rarefaction wave moves 
toward the center of the cylindrical geometry as a 
consequence the radial size of the emission zone de- 
creases with time. Therefore, the size of the emis- 
sion zone is larger at early times and smaller at late 
time. The high (M) regions are dominated by the 
early partonic phase where the collective flow has 
not been developed fully i.e. the ratio of collective 
to thermal energy is small hence show larger R S id e 
for the source. In contrast, the lepton pairs with 
M ~ TO p are emitted from the late hadronic phase 
where the size of the emission zone and hence the 
Rside is smaller due to larger collective flow. The 
ratio of collective to thermal energy for such cases 
is quite large, which is reflected as a dip in the vari- 
ation of Rside with (M) around the p-mass region 
(Fig. [3J upper panel) . Thus the variation of R S ide 
with M can be used as an efficient tool to measure 
the collectivity in various phases of the matter. The 
dip in R S ide at (M) ~ m p is due to the contribution 
dominantly from the hadronic phase. The dip, in 
fact vanishes if the contributions from p and u> are 
switched off (circle in Fig. [3]). We observe that by 
keeping the p and u contributions and setting radial 
velocity, v r = 0, the dip in R S ide vanishes, confirm- 
ing the fact that the dip is caused by the radial flow 
of the hadronic matter. Therefore, the value of R S ide 
at (M) ~ m p may be used to estimate the average 
v r in the hadronic phase. 

The Rout probes both the transverse dimension 
and the duration of emission. Therefore, unlike 
Rside, it does not remain constant even in the ab- 
sence of radial flow. The large M regions are pop- 
ulated by lepton pairs from early partonic phase 
where the effect of flow is small and the duration 
of emission is also small - resulting in smaller values 
of Rout- For lepton pair from M ~ m p the flow is 
large which could have resulted in a dip as in Rside in 
this M region. However, R ou t probes the duration 
of emission too which is large for hadronic phase. 
The larger duration compensates the reduction of 
Rout due to flow in the hadronic phase resulting in 
a bump in R ou t in this region of M (Fig. |3J lower 
panel). Both R S ide and R ou t approach QGP values 
for (M) ~ 2.5 GeV implying dominant contributions 
from the partonic phase in this region of M. 

It is clear from the results displayed in Fig. [5] 
that the high M(> rn ) and low M(< 0.5 GeV) 
regions are dominated by the radiation from early 
QGP phase when collective motion is not fully de- 
veloped - as a result the size of the emission zones 
for these M regions are large. However, lepton pairs 
for M around p mass dominantly originate from the 
late hadronic phase when the collectivity is large and 
consequently the size of the homogeneous emission 
domain is small. This results in non-monotonic vari- 
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FIG. 4: (color online) The ratio Rout /Rside and the dif- 
ference a/ R ut ~ Rside as a function of (M). 



ation of R S ide with M as discussed above. Therefore, 
such a non-monotonic behaviour of the HBT radii 
will signal the presence two different phases dur- 
ing the evolution - indicating the creation of QGP 
which will inevitably reverts to hadrons because of 
the cooling due to expansion. 

The quantities R ou t and Rside are proportional to 
the average size of the system [2J|. However, in the 
ratio, Rout/ Rside some of the uncertainties associ- 
ated with the space time evolution get canceled out. 
The quantity, Rout I Rside gives the duration of par- 
ticle emission |25l - [27| for various domains of M. 

Figure @] shows the variation of the ratio, 
Rout/ Rside and the difference, a/ R 2 out - R 2 slde as a 
function of (M) for Au+Au collisions at ^snn = 
200 GeV. Both show a non-monotonic dependence 
on (M). The smaller values of both the quantities, 
particularly at high mass region, reflect the contri- 
butions from the early partonic phase of the system. 
The peak around p-meson mass reflects dominance 
of the contribution from the late hadronic phase as 
discussed before. 

Now we discuss below two experimental challenges 
in such studies. We quote some numbers from the 
PHENIX measurements, keeping in mind that the 
situation will further improved by increasing the lu- 
minosity as well as collision energy. The number of 
events can be computed from the luminosity (£) , the 
pp in-elastic cross-section (cr) and the run time (T) 
of the machine as, 

N CV ent = Cx a xT (6) 
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For RHIC, the luminosity, C is of the order of 
50 x 10 27 cm _2 sec _1 and er=40 mb. If RHIC runs 
for 12 weeks then the number of events, N cvent = 
1.45 x 10 10 . For the M, 810 < M(MeV)< 990, the 
differential number (dN /2TTpxdpTdy) measured by 
PHENIX collaboration in Au+Au collisions [l(| at 
-y/SNN = 200 GeV is give by: 



dN 



2npTdpTdy 



\y = 0= ^± x 1.29 x 10- 7 (7) 



in the pr bin 1-1.8 GeV. Therefore, the (differential) 
number of pairs in the above range of px and M is 



1.45 x 10 



10 



x 1.29 x 10" 



N„ 



x 1870. 



Similarly for M, 500 < M(MeV)< 750 2 , the mea- 
sured value of the above quantity is: 



dN 
2irpTdpTdy 



N, 



\v=o~- 



part 



x 2.235 x 10 



-7 



(8) 



for the pt bin 1.4-1.8 GeV. This indicates that the 
(differential) number of pairs in this kinematic do- 



main is 



1.45 x 10 



10 



x 2.235 x 10~ 7 



x 3240. 



For — 10% centrality the number of participants 
for Au+Au collisions at ^/snw = 200 GeV is about 
330. The number of lepton pairs in the p? range 1.4- 
1.5 GeV is - 5.3 x 10 5 for the M window 0.5-0.75 
GeV. For 12 weeks of run time the number of events 
estimated with the current RHIC luminosity is ~ 
1.45 x 10 10 . Then the number of pairs produced per 



event is ~ 3 x 10 5 in the kinematic range mentioned 
above. The probability to have 2 pairs of dileptons is 
- 10~ 9 . Therefore, roughly 10 9 events are required 
to make the HBT interferometry with lepton pairs 
possible. 

It is expected that further increase in luminos- 
ity at RHIC by a factor 2 beyond the year 2012 to 
about 10 29 cm~ 2 s~* may be a motivating factor for 
such measurements. The increase in lepton pair pro- 
duction at higher collision energy at Large Hadron 
Collider (LHC) may also provide a reason to pursue 
such measurements. 

The possibility of dilution of signal due to addi- 
tion of random pairs, which one may encounter in 
the analysis of experimental data is discussed be- 
low. We have added some " mixture" to the dilepton 
source with exponential energy distribution, i.e. we 
have replaced u by u> + 8lo where 5uj has exponential 
energy (of the pair) dependence and weight factor is 
as large as that of ui itself. Then we find that the re- 
sulting change in the HBT radii is negligibly small. 
This can be understood from the expression for C 2 
(Eq.EQ: 



Co 



+ 



where 



^ J d 4 xiUj(xi,K)cos(ai) J d 4 X2^(x2,K) cos(a 2 ) 

/ d 4 xui(x, k\) J d 4 xuj{x,k2) 
J d 4 X\Uj(xi, K) sin(ai) J d 4 X2w(x2, K ) sin(a 2 )] 
J d 4 xuj(x,ki) J d 4 xui{x,k2) 



J 



q:i = tiMit cosh(?/i - t]i) — rikiT cos(#i — ipi) — t i M2T cosh(y 2 — Vi) + ?"1&2T cos(#i — ^2) 
012 = T2M2T cosh(y 2 - 772) - r 2 fc 2 T cos(# 2 — -02) — t 2 A/it cosh(yi — ?y 2 ) + r 2 fciT cos(# 2 — 
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where Xi — (TiCoshrji,riTCOs6i^riTCOs9i,TiSinhrii). 

It is clear that the expression for C 2 contains 
quadratic power of the source function both in the 
numerator and denominator. Therefore, changes in 
the source function will lead to some sort of partial 
cancellation (complete cancellation is not possible 
because the source function appears in the numer- 
ator and the denominator inside the integral with 
different dependent variables, K or fc^'s). 

Now at this point some comments on the change 
of the spectral function of light vector mesons in the 
medium are in order here. The modification of the 
spectral function of light vector mesons at non-zero 
temperatures and densities is a field of high con- 
temporary interest both from the experimental and 
theoretical (see 6, 8, 28] for reviews) point of views. 
The enhanced production of lepton pairs at the low 



M domain in In+In collisions at 158 A GeV beam 
energy measured by NA60 collaboration [9[ indicates 
substantial broadening of the p spectral function. In 
the present work we artificially modify the spectral 
function of the p by increasing its width by a fac- 
tor of two and keeping the pole mass at its vacuum 
value and then evaluate C 2 with the modified spec- 
tral function to explicitly check the change in the 
magnitudes of the HBT radii. We observe that an 
increase in the width of p by a factor of two leads 
to the reduction of R S ide by about 10% for invariant 
mass windows below p-peak. This reduction is due 
to the contributions from the decays of broader p 
which undergoes larger flow compared to the QGP 
phase. The nature of the variation of R S ide with (M) 
remains unaltered. However, a huge broadening of p 
for which the contribution from the hadronic phase 



G 



becomes overwhelmingly large compare to the QGP 
phase may alter the nature of the variation of R S ide 
in the low M domain. 



IV. SUMMARY 

In summary, the dilepton pair correlation func- 
tions has been evaluated for Au+Au collisions at 
RHIC energy. The additional kinematic variable, M 
for dilepton pairs make it a more useful tool for char- 
acterizing the different phases of the matter formed 
in heavy ion collisions compared to the HBT inter- 
fcromctry with direct photons. The HBT radii ex- 
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